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ABSTRACT 

We present a sample of 167 type 2 AGNs with double-peaked [Oiii] AA4959,5007 narrow emission lines, 
selected from the Seventh Data Release of the Sloan Digital Sky Survey. The double-peaked profiles can be 
well modeled by two velocity components, blueshifted and redshifted from the systemic velocity. Half of 
these objects have a more prominent redshifted component. In cases where the H/3 emission line is strong, it 
also shows two velocity components whose line-of-sight (LOS) velocity offsets are consistent with those of 
[Oiii]. The relative LOS velocity offset between the two components is typically a few hundred kms~', larger 
by a factor of ~ 1 .5 than the line full width at half maximum of each component. The offset correlates with 
the host stellar velocity dispersion cr,. The host galaxies of this sample show systematically larger c*, stellar 
masses, and concentrations, and older luminosity-weighted mean stellar ages than a regular type 2 AGN sample 
matched in redshift, [Oiii] A5007 equivalent width and luminosity; they show no significant difference in radio 
properties. These double-peaked features could be due to narrow-line region kinematics, or binary black holes. 
The statistical properties do not show strong preference for or against either scenario, and spatially resolved 
optical imaging, spectroscopy, radio or X-ray followup are needed to draw firm conclusions. 
Subject headings: black hole physics - galaxies: active - cosmology: observations - quasars: general - surveys 



1. INTRODUCTION 

Binary^ supermassive black holes (SMBHs) are pos- 
sible o utcomes of the hierarchical mergers of galaxies 
(e.g., Bege lman etall Il980t iMilosavlievic & Merrittj 120011: 
|Yu 2002). In one of the leading hypotheses, major merg- 
ers between galaxies are responsible f or triggering nu - 
clear starbursts and quasar activity (e.g., iHernquisl Il989h . 
Despite the success of the merger scenario in explain- 
ing much of the observed phenomenology of AGN statis- 
tics (e.g., Kauffmann & Haehnelt 2000; Volonteri et al. 2003; 
IWvithe & Loeh 2003; H opkins et all I2008t iShen 2009) and 
the cores of massive elliptical galaxies (e.g.. lRiber et al.il99^ 
iKormendy & Bender 2009), direct observational evidence for 
binary SMBHs is surprisingly scarce. 

The fraction of binary quasars at separations of tens to hun- 
dreds of kpc (h alo) scal es, is 5, 0.1% at 1 5. z 5. 5 (e.g ., 
[Hennawi et al.l2006; Mve rs et aU2008l:lHennawi et al.l2009l) . 
On kpc (galactic) scales there are only a handful of un- 
ambiguous low-redshift cases known, in which both active 
SMBHs are detected in X-rays (NGC 6240 and Mrk 463, 
iKomossa et al."2003'; 'Bianchi et al.'l 200l) . or in the radio (3C 
75, Owen et al. 1985). On sub-kpc scales, only one case is 
known (0402-H379. lRodriguez et al. 2006) of a pair of BHs de- 
tected by VLBI with a projected separation of ^ 7 pc. There 
are several possible reasons that binary SMBHs are so rare 
compared to the expectations from the merger scenario: they 
spend most of their time at separations far below kpc scales 
and hence are difficult to resolve spatially; one or both of the 
SMBHs are heavily obscured; both BHs are rarely active at 
the same time, and dynamical differences between a single 
BH and a binary BH are not readily discernable at cosmolog- 
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ical distances; 

Recently, IComerford et al.l (12009 a|) conducted a systematic 
survey of typ e 2 AGNs in the DEEP2 galaxy sample at z ~ 0.6 
(iDavis et all 12003 1). They argue that 37 of the 107 AGNs in 
their sample are inspiralling binary SMBHs, based on signif- 
icant velocity offsets between the narrow-line region (NLR) 
emission lines [O iii] AA4959,5007 and the systemic redshift 
measured from stellar absorption fe atures. In particul ar, two 
of these objects (one reported by iGerke etaP l2007h show 
double-peaked [O iii] lines with velocity offsets of a few hun- 
dred kms"' and spatial offsets of several kpc; they interpreted 
these as binary SMBHs when both BHs are active. Such 
double-peaked [O iii] emission line features have also been 
seen i n a few optically-selected AGNs (e.g.. | Heckman et al.l 
[19811 IZakamska et al.H200l IZhou et al.ll2004l) . One advan- 
tage of using spectroscopy to find binary SMBHs is that it is 
not limited by spatial resolution, as long as the separation is 
not so small that the NLRs are no longer distinct. On the other 
hand, such spectral features may be due to biconical outflows 
or disk rotation around a single SMBH (e.g. , lAxon et al.l 19981; 
IVeilleux et alll200Tt ICrenshaw et al.ll2009l) . Here we focus 
on double-peaked narrow line objects in which the associ- 
ated line emitting gas is on scales of ^^100 pc to several kpc. 
Therefore in the binary SMBH scenario, our double-peaked 
sample would include objects whose two SMBHs along with 
their own NLR gas are still co-rotating in the galactic po- 
tential, well befo re forming a pc-scale binary under dynam- 
ical friction (e.g.. lMilosavljevic & Merrit3l200ll) . For double- 
peaked broad line objects whi ch may involve outflows or 
dis ks on 5,p c scales, see, e.g., lEracleous & HalpernI (119941) 
and lStrateva et al. (2003,) . 

We have carri ed out a systematic search in the SDSS 
dYork et al.l 120001 ) spectroscopic database for emission line 
AGNs with double-peaked features or peculiar [Oiii] line 
shifts relative to the systemic redshift as measured from stel- 
lar absorption features. Here we report our initial results on a 
sample of 167 type 2 AGNs with double-peaked [Oiii] emis- 
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TABLE 1 
The sample 



SDSS Designation Plate Fiber MJD z ct* FWHM[o iii], i FWHM[oiii].2 ^'[om].i V[oin].2 Vhsj Vh/3.2 



J000249.07 +004504.8 388 345 51793 0.0868 243 271 219 -361 168 -282 197 

J00091 1.58 -003654.7 388 148 51793 0.0733 203 231 220 -198 134 -196 129 

J010750.48-005352.9 670 55 52520 0.5202 281 645 -546 37 -540 -113 

JOl 1659.59 -102539.1 660 213 52177 0.1503 181 304 203 -159 145 -81 189 



Note. — The full table is available in the electronic veision of the paper. The subscripts ''1" and "2"' denote blueshifted and redshifted components. K[o ni] and Vnf^ 
denote velocity offsets relative to the systemic redshift measured over [O III] A A4959,5007 and H/3. Velocities are in units of km s"' . The entry reads zero if the quantity is 
unmeasurable. The typical uncertainty in is ~ 15kms"' . For [O III], the typical statistical enors in the best-fit velocity offsets and FWHMs are ~ 5 km s"^ and ~ 10 km 
s * . For H/3, the statistical errors are generally a few times larger because H/3 is weaker than [O III]. 
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Fig. 1 . — Redshift distribution of the parent AGN and double-peaked sam- 
ples. 

sion lines^. We describe our sample selection in ^ the sta- 
tistical properties of the sample are presented in ^ and we 
discuss our results and conclude in ^ A cosmology with 
0,„ = 0.3, Ha = 0.7, and h = 0.7 is assumed throughout. 

2. DATA 

2.1. The Sample 

Our parent sample is the MPA-JHU SDSS DR7 galaxy 
sample^ drawn from t h e SD SS DR7 spectroscopic 
database (lAbazajia n et al.l 120091) : those include objects 
spectrally classified as galaxie s by the specBS pipeline 
(lAdelm an-McCarthv et al. 2008) or quasars that are targeted 
as galaxies ( Strauss et al. 2002; Eisenstein et al. 2001) and 
have redshifts z < 0.7. We adopt redshifts z and stellar 
velocity dispersions a, from the specBS pipeUne. When 
unavailable from the pipeline, we measure us i ng the 
direct-fitting algorithm described in ' Greene & Hd (120061) 
and |Ho et al. (2009); we have checked that this code gives 
results consistent with those of specBS. We have checked 
the accuracy of the redshift by refitting the stellar continuum 

^ After we submitted the paper, two similar statistical studies by 
ISmith et al] ^009) and Wans et al. (2009) appeared on arXiv. Smith et al. 
(20090 focus on broad line AGNs with double-peaked [O III] lines selected 
from SDSS quasar spectra and find a similar fraction (~1%) of double- 
peaked objects that we found. Because of pipeline misclassification, their 
double-peaked sample also con tains narrow line objects, some of which over- 
lap with our catalog. Wang et ajj j2009l) studied a sample of 87 double- 
peaked narrow line AGNs of which 42 objects overlap with our sample. 

' http://www.mpa-garching.mpg.de/SDSS/ 



with galaxy templates (iLiu et al.ll2009t) . as well as comparing 
multiple epoch spectra of the same objects. The quoted 
redshift faithfully traces the stellar absorption features (which 
we take as the systemic redshift), with typical statistical 
errors of a few kms"' and systematic errors of ~ 10 km 
s"'. Additional spectral and photometric properties such 
as emission-line fluxes (determined from equivalent widths 
and continuum fluxes), spectral indices, stellar masses, and 
concentration indices are taken from the MPA-JHU SDSS 
DR7 data product. 

We select our AGN sample from this parent sample 
by the following criteria: 1) the rest-frame wavelength 
ranges [4700,5100] A and [4982,5035] A centered on the 
[O III] A5007 line have median signal-to-noise ratio (S/N) > 5 
pixel"' and bad pixel fraction < 30%; 2) the [Oiii]A5007 
line is detected at > 5a and has a rest-frame equivalent width 
(EW) > 4 A; 3) the Hne flux ratio [Oiii] A5007/H/3 > 3 if 
z > 0.33, or the diagnostic Hne ratios [Oiii] A5007/H/3 and 
[Nil] A6584/Ha lie above t he theoretical upper U mits for star- 
formation excitation from 'Ke wley et al.l (1200 I h on the BPT 
diagram (Baldwin et al. 1981) if z < 0.33. This procedure 
yields ^ 14,300 type 2 AGNs. We supplement this sample 
with ^ 400 type 2 quasars from Reves et al. (2008) which are 
not included in the MPA-JHU data products. These additional 
type 2 quasars extend to somewhat higher redshifts (z < 0.83). 
Our final AGN sample includes 14,756 objects with high S/N 
and high spectral quality around the [O iii] lines, suitable for 
the analysis that follows. The redshift distribution of this par- 
ent AGN sample is shown in Figure [T] 

In many AGNs and starbursts, the narrow forbidden emis- 
sion lines [O iii] AA4959,5007 are known to have peculiarities 
such as extended blue wings, velocity offsets from systemic 
redsh ifts, and complex line profiles (e.g ., Heckman et aD 
TMT; WhittleT985a:' Zakamska et al.l2003tlZhou et aLIBOO^ 
Komossa et al. 2008). We here focus on a more dramatic 
subset - those with double-peaked [Oiii] emission lines. A 
complete analysis of the nature and statistical properties of 
other peculiar [Oiii] line characteristics will be presented 
elsewhere. 

The identification of such double-peaked objects requires 
that both [O iii] A4959 and [O iii] A5007 are better fit with two 
components rather than a single component. As the initial 
screening process we visually inspected the spectra of all the 
14,756 AGNs and identified 167 objects with unambiguous 
double-peaked [Oiii] lines. We only include objects which 
have well-detected double peaks in both [O iii] A4959 and 
[Oiii] A5007 with similar profiles; we do not include those 
with complex line profiles such as lumpy, winged, or multi- 
component features. We plot the redshift distribution of the 
double-peaked sample in Figure[T] and the SDSS images and 
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Fig. 2. — Example double-peaked AGNs. Shown here are SDSS gri color-composite images and host-subtracted spectra (data points in red and smoothed 
curve in black) along with our best fits (model in green and individual components in magenta) for the H/3-[0 III] region. The vertical lines are drawn at the 
systemic redshift of each galaxy as measured from stellar absorption features. 



spectra of three such examples in Figure |2l We Hst all the 
167 double-peaked AGNs and their line measurements in Ta- 
ble [T] We caution that there could be some double-peaked 
narrow-line objects in SDSS DRV missed by our selection if 
the pipeline got the redshifts wrong because of the peculiar 
line profiles, although this fraction is likely to be no larger 
than a few percent. 

We measured the line properties of these double-peaked ob- 
jects as follows. The galaxy continuum was subtracted using 
the best-fit template constructed fr om a linear combination of 
instantaneous starburst models of iBruzu al & Chariot (2003) 
with ten different ages, as described in iLiu et al., (.2009,) . In 
the 1 1 objects for which the stellar continuum is too weak 
for measuring cr* and for template fitting, a simple power- 
law model was used. The continuum-subtracted [O iii] region 
(AA4930-5040A) was fit by a pair of Lorentzian functions 
convolved with the measured instrumental resolution of the 
spectra (a ~ 65 km s~'). The redshift and line width for 
each velocity component of [O iii] A4959 and [O iii] A5007 



were forced to be the same. The flux ratio of [O iii] A5007 
to [O III] A4959 was allowed to vary (but we found it was al- 
ways close to 3). In cases where H/3 is measurable, we further 
fit a double-Lorentzian to the Hf3 region (AA4850-4880 A), 
where the positions of the centroids were allowed to vary, but 
the widths of the two H/? components were fixed to the best- 
fit values for the two [Oiii] components. We then repeated 
this fitting process with a double-Gaussian model, and took 
as the best fit the model which had the smaller reduced for 
each object. We visually inspected all the fits and verified that 
the line profiles are generally well reproduced by the model. 
For [Oiii], the typical statistical errors in the best-fit velocity 
difference between the two components Voff, line fluxes and 
FWHMs are 5 km s~' , ~ 5%, and ~ 10 km s~' , but we cau- 
tion that the actual uncertainties are likely to be larger For 
H/3, the statistical errors are generally a few times larger be- 
cause H(3 is weaker than [Oiii]. 

2.2. Selection Completeness 
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Fig. 3. — Selection completeness as a function of median S/N, v^ff, FWHM, total [OIII] A5007 EW and the flux ratio between the two components, as 
determined from our Monte Carlo simulations. In the middle panel, black and blue hues are for the two components respectively. 
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Fig. 4. — Narrow-line region physical properties of the double-peaked components. Objects with stronger redshifted components are plotted as red open 
circles whereas those with stronger blueshifted components are plotted as blue plus signs. The Spearman correlation coefficient and the probability of null 
correlation are labelled in each plot for the whole (black), the stronger-blueshifted (blue), and the stronger-redshifted (red) samples, respectively, (a). Ionization 
parameter (indicated by the flux ratio between [O III] A5007 and H/3) of the blueshifted component versus that of the redshifted component, (b). Line width 
ratio of the double-peaked components versus Balmer decrement (as a measure of reddening), (c). Electron density (inferred from the flux ratio between the 
[S II] AA6717.6731 doublet) of the blueshifted component versus that of the redshifted component. 



Our selection of double-peaked objects was done by eye 
and is by no means objective. To gain some sense of our com- 
pleteness we performed Monte Carlo simulations. We gener- 
ated mock spectra of double-peaked objects in the H/3-[Oiii] 
region with random distributions of continuum S/N, EW, ve- 
locity offset, and FWHM of the two components of [Oiii]. 
We randomly assigned a Gaussian or Lorentzian profile to 
each line and broadened the spectrum using an instrumental 
Gaussian broadening of = 65 km s"' . In this way we gener- 
ated 2,000 mock spectra of double-peaked objects and mixed 
them with 18,000 mock spectra of single-peaked objects with 
random spectral properties. We then visually inspected the 
20,000 mock spectra and identified double-peaked objects in 
the same way as we did for the real sample. We show in Fig- 
ure[3]the selection completeness as a function of median S/N, 
Voff, FWHM, total [Oiii] A5007 EW, and flux (EW) ratio be- 
tween the two components /ew- In general the completeness 
increases with increasing Voff, [Oiii] A5007 EW, and /ew, and 
decreases with increasing FWHM, as expected from the way 
our double-peaked objects are identified. As a result there is a 
selection bias against objects with small Voff and large FWHM 
that we discuss later. In addition we are very incomplete for 
objects with Voff < 200kms~' partly due to the spectral res- 
olution limit*'. On the other hand, we find that the false de- 

For instance, the object SDSS J100043.13-H020637.2, with v„ff ~ 
150kms"' reported in lComerford et all <2009bl) is not included in our sample 
based on the SDSS spectrum alone. 



tection rate is very low (< 2% of the double-peaked sample) 
- probably an advantage of visual inspection over automated 
algorithms. We note that the completeness estimated this way 
is not the absolute completeness, which depends on the actual 
underlying distributions of all relevant properties. Neverthe- 
less, these completeness estimates are useful to correct for se- 
lection biases. For example, when other properties are fixed, 
double-peaked objects with larger EWs are easier to identify, 
by a relative amount that can be determined from our simula- 
tion results. 

After accounting for the dependence of selection complete- 
ness on [O III] EW, we find tentative evidence that the fraction 
of double-peaked objects in the underlying AGN sample de- 
pends on [Oiii] A5007 luminosity (L[Oiii]): it increases by a 
factor of ^ 2 from L[Oiii] ^ lO^ ^L© to L[Oiii] ~ lO^'^L©, and 
has no clear dependence on redshift in the range probed. In 
most cases the two velocity components are redshifted and 
blueshifted with respect to the systemic redshift (except that in 
a few cases, one component is apparently coincident with the 
systemic velocity within the errors), and about half of them 
have more prominent redshifted components - a subset which 
may be more unusual since objects with outflows may tend to 
have more prominent blueshifted components (see below). In 
cases where H/3 is measurable, it also shows a double-peaked 
feature with a similar velocity offset as that of [Oiii]. Here- 
after we use Voff measured from [O iii] in our analysis as it is 
more robust than that from H/3. As illustrated in Figure |4j a), 
both velocity components have excitation diagnostic line ra- 
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Fig. 5. — Host galaxy properties (stellar velocity dispersion ct*, stellar mass M*, spectral indices D„(4000) and H5a indicating the luminosity-weighted mean 
stellar age and the post-starburst fraction, and r-band concentration index Rqo.r/^so.r) of the double-peaked sample (red) compared to the control sample (black). 
The control sample is randomly drawn from the parent AGN sample to have identical redshift and [O III] A5007 EW distributions as the double-peaked sample. 
For each distribution the median is marked with a vertical line. 



tios H/3/[Oiii] characteristic of AGNs^ The FWHM of each 
velocity component is typically smaller than Voft by a factor 
of - L5. 

3. STATISTICAL PROPERTIES 

3.L Host Galaxy Properties 

Are there any characteristics of the host galaxies of the 
double-peaked objects that are different from ordinary AGNs? 
The double-peake d sa mple is not selected uniformly as we 
demonstrated in ^2.21 To account for selection biases, we 
construct a control sample from the underlying AGN sam- 
ple with identical redshift, [O iii] A5007 EW and luminos- 
ity distributions as the double-peaked sample. We compare 
their host-galaxy properties, including stellar mass and ve- 
locity dispersion, SDSS colors, stellar population parame- 
ters D„(4000) and HSa (indicating the luminosity-weighted 
mean stellar age and the post-starburs t fraction in the past 
— 0.1-1 Gyr, respectively; e.g.. .Kauffmann et alj 120031) . in- 
trinsic extinction (estimated by the Balmer decrement), effec- 
tive radius, concentration (characterized by, Rgoj/R^o.r, the 
ratio of radii that contain 90% and 50% of the r-band Pet- 
rosian flux), and inclination (estimated by the ratio between 
the major and minor isophotal axis). As shown in Figure |5] 
the double-peaked sample has systematically larger stellar ve- 
locity dispersions and masses, older mean stellar ages (and 
redder colors), smaller fractions of post-starburst populations 
in the past 0.1-1 Gyr, and higher concentrations than the con- 
trol sample. The difference in the centroids of the two dis- 
tributions are roughly half of the width of the distributions. 
Kolmogorov-Smirnov (KS) tests show that the probabilities 
that the two samples are drawn from the same distribution are 
Pks < 1Q~^. All other properties of the host galaxies we exam- 
ined were indistinguishable with KS probabilities Pks > 10"'. 
We also match t he double-peaked and the control sa mples 
with the FIRST (iBecker et all [1991 IWhite et all [1997 ) and 
ROSAT (iVoges et all II 999L 120001)" surveys and find that they 
contain very similar fractions of matches with both surveys 
(- 35% with FIRST and - 2% with ROSAT). The subset that 
has a more prominent redshifted component is indistinguish- 
able from the rest of the double-peaked sample in terms of the 
properties studied above. 

3.2. Correlations between Line Properties 

There are several apparent correlations among the dy- 
namical properties and the NLR physical conditions of the 

' For objects with z < 0.3, the [N II] and Ho lines are visible; they appear 
double-peaked as well, although the different components often are blended 
with one another. 



double-peaked sample. Figures Ua) and |4jc) illustrate that 
the two velocity components have similar ionization param- 
eters (as indicated by the diagnostic line ratio H/3/[Oiii]) 
and electron densities (measured from the ratio of the dou- 
blet [Sii] AA6717,6731 for 107 objects in the sample which 
have low enough redshifts to have [S ii] coverage in the SDSS 
spectra and whose inferred electron densities of both compo- 
nents fall in the range of [10, 10"^] cm"^). Figures |6ja)-(d) 
show that Voff is correlated with both ct* and the line FWHMs 
of the two components, and the line FWHMs are correlated 
with each other. In all these four relations, the correlation for 
objects with more prominent redshifted components appears 
somewhat stronger than that for the rest. There is some selec- 
tion bias at the low Voff versus large FWHM corner, where the 
double-peaked feature is difficult to identify. However, our 
Monte-Carlo simulations show no selecti on b ias against ob- 
jects with high Voff and small FWHMs ( ^2. 2b . In particular, 
the correlation between Voff and ci* suggests that the double- 
peaked line-emitting regions are on galactic scales where the 
galactic potential dominates the bulk kinematics. Further- 
more, Figures |6je) and|6jf) show that both the [O iii] line flux 
ratio and the line width ratio of the two velocity components 
are anti-correlated with the ratio of their L OS velocity offset s 
relative to the systemic redshift (see also IWang et al.ll2009h . 
which most likely results from momentum conservation. 

4. DISCUSSION 

4.1. NLR Kinematics 

The double-peaked feature may result from particular NLR 
geometries such as biconical outflows** or rotating disks on 
kpc scales. In these scenarios, there is only one SMBH 
and the observed blueshifted and redshifted emission-line 
peaks arise from NLR gas moving towards and away from 
us. Certain nearby Seyfert galaxies known to have biconi- 
cal outflows show double-peaked emission when spectra are 
taken over the whole galaxy. Such p rototypical cases in- 
clude NGC 1068 (e.g..lAxon et al.ll998h and NGC 3079 (e.g., 
Duric & Seaquisi [T988l) . It is very likely that there are some 
such sources in our sample. Several objects in our sample 
show clear galactic disks in their SDSS images. Here the cor- 
relation between line width and velocity shift naturally arises, 
as the bulge gravitational potential well influences the motion 
of NLR gas. The correlation between the diagnostic line ratio 
[O iii]/H/3 of the two components shown in Figure |4ja) may 

In principle, a NLR with inflowing gas and dust can also produce 
the observed features, but this scenario seems less favored considering the 
difficulty of removing high angular momentum given the large v^ff (e.g., 
iHeckman et al.iI98Ii) . 
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Fig. 6. — Dynamical properties of the double-peaked components. Objects with stronger redshifted components are plotted as red open circles whereas those 
with stronger blueshifted components are plotted as blue plus signs. The Spearman correlation coefficient and the probability of null correlation are labelled in 
each plot for the whole (black), the stronger-blueshifted (blue), and the stronger-redshifted (red) samples, respectively, (a)-(c). stellar velocity dispersion cr*, 
FWHM of the blueshifted component, and FWHM of the redshifted component versus velocity offset between the double-peaked components, (d). FWHM 
of the blueshifted component versus that of the redshifted component, (e) Velocity-offset ratio of the double-peaked components versus the flux ratio, (f). 
Velocity-offset ratio of the double-peaked components versus the FWHM ratio. 



be explained if the NLR gas of both velocity components is 
ionized by a single SMBH. 

Can all the objects in our double-peaked sample be ascribed 
to NLR geometry? Outflows and inflows of NLR gas have 
long been invoked to explain several commonly observed fea- 
tures in NLR high ionization lines (such as [Oiii]), including 
extended wings (usually blueshifted from the systemic red- 
shift), line asymmetries, and broader line widths than those of 
low ionization Unes (e.g. , Heckman et al. 1981; Penston et al. 
119841; iRodn'guez-Ard ilaet al. 2006). In the cases o f NGC 
1068 ( e.g..lAxon et al . 1998) an d Mrk 78 (e.g., He ckman et alj 
[19811 lPedlaretal.1 11989; Whi ttle & Wil son 2004J) known to 
have biconical outflows, the velocity splitting relative to the 
systemic redshift is larger on the blueshifted component, 
which is also brighter and has a larger width than the red- 
shifted component. The st atistical signifi cance of these trends 
is unclear, however (e.g.. lWTiittle et al] |1988). Unlike these, 
our double-peaked sample has comparable line widths in the 
blueshifted and redshifted components as shown in Figure 
|6fd); half of our double-peaked objects have a more promi- 
nent redshifted component than the blueshifted component. 
In addition, our sample has consistent Vott in H/3 and in [O iii] 
whenever H/3 is measurable, with no significant ionization 
stratification as expected in NLR outflows. These possible 
discrepancies may arise because our sample has appreciably 
larger [O iii] luminosity than these previous low redshift com- 
parison samples, and these properties may be a function of lu- 



minosity or redshift or both; alternatively, they may suggest 
that not all objects in our double-peaked sample exhibit NLR 
outflows. 

Further possible clues come from host galaxy and line prop- 
erties. The FIRST detected fraction for the double-peaked 
objects is comparable to that of the control sample. The host- 
galaxy inclination distribution shows no significant difference 
either, whereas our selection of a large LOS v^ff for the double 
components would bias towards either edge-on or face-on in- 
clinations for the rotating disk or bipolar outflow populations, 
unless the outflow orientation is not closely linked to the an- 
gular momentum axis of the host galaxies. Radio outflows in 
nearby Seyfert galaxies and low-ionization AGNs indeed ap- 
pear to orie nt randomly with respect to the host gala xy axes 
(e.g., Ulves tad & Wil son' 1984': iGalUmore e"ta r2006). In ad- 
dition, while a correlation bet ween Vpff and line FWHM is ob- 
served in NLR outflows (e.g.. lKomossa e t al. 2008), the ratio 
of velocity shift to line width is of order unity, unlike the ratio 
of ^ 1.5 seen in our double-peaked sample. Finally, as shown 
in Figure Etb), there is no obvious correlation in our sample 
between the Balmer decrement and the line-width ratio be- 
tween the two components of the double-peaked feature (as a 
proxy for asymmetry). Such a correlation might exist if the 
blueshifted and redshifted NLR gas components are intrin- 
sically symmetri c and the appare nt asymmetry results from 
extinction. Hec kman et al.l (Il981h observed a correlation be- 
tween the asymmetry of a single line profile (which, how- 
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ever we caution, is different from our asymmetry indicator 
- the line width ratio of the double-peaked components) and 
Balmer decrement in a sample of 36 nearby Seyferts thought 
to co ntain outflows, although later studies on nearby Seyferts 
(e.g.. lWhittldll985"bl) found no strong correlations. 

4.2. Merging SMBH Pairs 

Another explanation for the double-peaked features is 
merging SMBH pairs. In this scenario, both BHs, along 
with their own NLR gas (with scales of order hundreds of 
pc), are co-rotating in the galactic potential on > kpc scales, 
well before dynamical friction causes their orbit to decay 
and form a gravitationally bounded compact binary (e.g., 
pilosavljevic & Merritt 2001). The Voff-cr* and Voff-FWHM 
correlations in this scenario would naturally arise from dy- 
namics. The comparable numbers of objects with more 
prominent redshifted and blueshifted components and the null 
correlation between the Balmer decrement and the double- 
peak asymmetry may both be accommodated if the two emis- 
sion components are associated with two SMBHs. The simi- 
lar distribution of host-galaxy inclination angles to the control 
sample could be explained, given that there is no preference 
for the rotation axis of galaxies in the merging SMBH-pair 
scenario. 

On the other hand, the correlation between the diagnostic 
line ratio [Oiii]/H/3 of the two components shown in Fig- 
ure Ua) may be difficult to explain if the NLR gas of each 
component is primarily ionized by its own BH. In addition, 
the higher galaxy concentration for objects with double peaks 
(Figure |5]l seems inconsistent with these objects being in the 
early stage of mergers, although we may be biased to more 
compact objects by the necessity to get both merging objects 
within the 3-arcsec fiber aperture. The older mean stellar ages 
and the smaller post-starburst fraction in the past 0.1-1 Gyr 
also seem counterintuitive, if there is excess star formation in 
mergers. However, these diffe rences might also resul t from 
the mass-age correlation (e.g.. iKauffmann et al.ll2003h given 
that the median cr* of the double-peaked sample is larger than 
that of the control sample (Figure |5]l, or if the associated star 
formation is on a different time scale from what is probed by 
the D„(4000) and HSa indices adopted here, or it is suppressed 
by AGN activity. 

In view of the above arguments, we conclude that our 
double-peaked sample may contain a population of NLR out- 
flows or rotating disks, as well as objects that are merging 
SMBHs; The statistical properties of the sample do not show 
any strong preference for or against either scenario. In Fig- 
ure |2] we show several intriguing cases. The first example 
(SDSSJ114642.5+511029.9) is a merging system of which 
the off-centered galaxy is partially covered by the 3-arcsec di- 
ameter fiber. We have checked tha t the two galax ies are both 
detected in the band by 2MASS (iSkrutskie et al.ii2006) . and 
have found evidence that both are AGNs using spatially re- 
solved longslit spectroscopy, the results of which will be pre- 
sented elsewhere. Only 13 of the 167 objects in our double- 
peaked sample show possible evidence for two cores in the 
SDSS images within the 3-arcsec diameter of the spectro- 



scopic fibers. The second object (SDSSJ122709.8+124854.5) 
does not have a resolved double core in its SDSS image, and 
its emission lines have well separated peaks, symmetrically 
shifted from the systemic velocity, and the redshifted com- 
ponent is more prominent. High resolution HST imaging 
may help rev eal potential double core s unresolved in SDSS 
images (e.g., IComerford et alj l2009bl) . For example, one 
of our double-peaked objects (SDSSJ130128. 8-005804.3, at 
z = 0.2455) shows no double cores in its SDSS image, but 
'Zaka mska et al.l (12006 ) used HST imaging to show that it con- 
sists of two galaxies with a projected separation of 1.3 arcsec 
(corresponding to 5.0 kpc) which well fits into the SDSS 3- 
arcsec diameter fiber. We will need spatially resolved spec- 
troscopy to determine whether both of the galaxies are AGN. 
The last object in Figure|2]has a strong component coincident 
with the systemic redshift and a weaker blueshifted compo- 
nent and an even weaker redshifted component (not fitted), 
which could arise from a classic NLR region plus bipolar out- 
flows. Spatially-resolved optical imaging, spectroscopy, ra- 
dio and/or X-ray followup are still needed to help draw firm 
conclusions on the nature of these double-peaked narrow line 
objects. 
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